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We consider a harmonically driven acoustic medium in the form of a (finite length) highly nonlinear 
granular crystal with an amplitude and frequency dependent boundary drive. Remarkably, despite 
the absence of a linear spectrum in the system, we identify resonant periodic propagation whereby 
the crystal responds at integer multiples of the drive period, and observe that they can lead to 
local maxima of transmitted force at its fixed boundary. In addition, we identify and discuss 
minima of the transmitted force (“antiresonances”) in between these resonances. Representative 
one-parameter complex bifurcation diagrams involve period doublings, Neimark-Sacker bifurcations 
as well as multiple isolas (e.g. of period-3, -4 or -5 solutions entrained by the forcing). We combine 
them in a more detailed, two-parameter bifurcation diagram describing the stability of such responses 
to both frequency and amplitude variations of the drive. This picture supports an unprecedented 
example of a notion of a (purely) “nonlinear spectrum” in a system which allows no sound wave 
propagation (due to zero sound speed: the so-called sonic vacuum). We rationalize this behavior in 
terms of purely nonlinear building blocks: apparent traveling and standing nonlinear waves. 

PACS numbers: 45.70.-n 05.45.-a 46.40.Cd 


Introduction. In media that bear a crystalline lat¬ 
tice structure (perfect or imperfect), an understanding 
of the underlying linear spectrum is of paramount im¬ 
portance both for their linear wave transmission proper¬ 
ties [T] and for the emergence of anharmonicity-induced 
excitations such as discrete breathers [2]. However, in 
recent years there has been an increased interest in sys¬ 
tems that do not possess a linear spectrum. This stems 
to a considerable degree from the interest in discretizing 
prototypical nonlinear systems such as the Burgers equa¬ 
tion [3], but also from the intrinsic mathematical (and 
potentially physical) interest in exploring systems with 
purely nonlinear dispersion; a remarkable feature of the 
latter has been the identification of discrete analogues of 
compactons [1] and of compactly supported breathers [5] . 
A fundamental question that then emerges in the lat¬ 
ter context is whether an analogue of acoustic transmis¬ 
sion/spectral features exists in these highly (or purely) 
nonlinear systems. 

Importantly, the last two decades have also brought 
forth a remarkably accessible (theoretically, numerically 
and experimentally) example of a (tunably) highly non¬ 
linear system. This is in the form of granular crys¬ 
tals, consisting of lattices of beads interacting (chiefly) 
elastically, by means of the so-called Hertzian contact 
forces F = ASP. Here A is an elasticity and geom¬ 
etry dependent prefactor, S is the compression-induced 
displacement and p the nonlinear force exponent, typi¬ 


cally equal to 3/2. These acoustic systems enjoy a re¬ 
markable degree of tunability from homogeneous (same 
types of beads) to heterogeneous or disordered, and im¬ 
portantly for our current purposes from essentially lin¬ 
ear to weakly or even highly or purely nonlinear EHD- 
The latter aspect is controlled by the amount of a pri¬ 
ori compression (dubbed precompression i5o) externally 
imposed on the beads. When this precompression is ab¬ 
sent, i.e.. So = 0, the Hertzian interaction of the particle 
chain manifests its essentially nonlinear character. It is 
in that case that the speed of sound c = (Api5g~^)^/^ 
in the medium vanishes, the so-called sonic vacuum sce¬ 
nario. In this setting, highly nonlinear solitary waves, 
pioneered by Nesterenko Hilo], have been identified and 
are the dominant energy carriers in the system m- 

Our principal question is thus formulated at the junc¬ 
tion of these two important themes. On the one hand, 
for linear systems, acoustic wave propagation and the 
linear “phonon” spectrum are crucial for a fundamental 
understanding of propagating vs. evanescent waves in 
the system. On the other hand, no such entities arise 
in the sonic vacuum. In the latter case, highly local¬ 
ized (in fact, doubly exponentially localized m) solitary 
waves constitute the prototypical excitation of the sys¬ 
tem. Gan then one seek spectral properties/spectral re¬ 
sponse in such highly nonlinear (finite lattice) settings 
and, if the answer is positive, how does such a response 
emerge? Surprisingly, in our system of choice (towards 
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exploring these ideas), namely the granular crystal, we 
computationally illustrate the existence of such a “non¬ 
linear spectrum” as it arises in response to externally 
driving the system at different frequencies. The prin¬ 
cipal mechanism for the emergence of this spectrum is 
the (nonlinearly induced) periodic and potentially quasi- 
periodic response of the system depending on the rela¬ 
tive phases between the applied harmonic excitation and 
waves (“pulses”) that appear to propagate in opposite 
directions, or to form standing wave patterns within the 
medium. These inherently nonlinear notions have the 
potential to be widely applicable to systems without a 
linear spectrum, well beyond the specifics of our particu¬ 
lar paradigm, and hence are of broad and diverse interest. 

Only a few works have discussed the responses of gran¬ 
ular crystals to harmonic excitation. In |13| nonlinear 
resonances in precompressed granular crystals were stud¬ 
ied, together with certain bifurcations related with these 
motions; however, in distinction from the system consid¬ 
ered in the present work, the systems in m possess a 
linear spectrum due to applied static precompression. In 
m the forced response of a two-bead granular crystal 
was considered and the pass and stop bands of this sys¬ 
tem were confirmed; given, however, the light damping of 
that system, no clear resonance peaks could be detected 
since the high-amplitude response became chaotic. To 
date the resonance and antiresonance structure of non¬ 
linear sonic vacua have not been systematically studied; 
this is our aim here. 

Setup. We consider a chain of beads interacting elasti¬ 
cally via Hertzian contacts that are mathematically de¬ 
scribed (upon a suitable rescaling of amplitudes and of 
time m) as: 

Xi = {x,_i - - (Xi - a;i+i)+^ 

-l-A [(ij_i - Xi)R[xi-x - Xi) 

- (xi - x^+i)H(xi - x^+i)], (1) 

with i = 1,..., N with N = 11). The Xi’s here represent 
the displacement of the particle centers from their equi¬ 
librium positions. The first two terms on the right hand 
side capture the bead-bead elastic interaction in our finite 
lattice, while the last two account for the linear dissipa¬ 
tive dynamics, following the proposal of m- H denotes 
the Heaviside function, while (•)+ = max(',0). These 
terms appear because of the tensionless character of the 
system, which implies there are no forces when beads are 
not in contact. The rescaling used to bring the system 
in the form 0 together with the numerical values of the 
system parameters are given in the Appendix. 

The boundaries of the system under consideration 
consist of (a) an actuator exciting the (zero-th) bead 
with prescribed harmonic displacement at the left end 
and (b) an immobile wall on the right end. The lat¬ 
ter is implemented as a fixed additional bead for which 
xn+i = xn +1 = 0, while the former is associated with 


xq = Asin(a;t), i.e., a periodic excitation of a given fre¬ 
quency uj with an amplitude A. The result of this har¬ 
monic drive in the presence of linear coupling between 
adjacent nodes of a (finite) lattice would be to excite the 
quantized (based on the boundary conditions) wavenum¬ 
bers allowable by the lattice, and, accordingly, the fre¬ 
quencies corresponding to these wavenumbers based on 
the linear dispersion relation. In the present case, how¬ 
ever, such a linear dispersion relation is completely absent 
(and so is a linear spectrum of frequencies) due to the es¬ 
sentially nonlinear character of the examined sonic vac¬ 
uum. Accordingly, the resulting dynamics are expected 
to be completely tunable with the amplitude of the forc¬ 
ing, whereas the existence and structure of resonances in 
this system (i.e., of amplitude amplification on the gran¬ 
ular chain at certain frequencies) is an interesting ques¬ 
tion. In the absence of linear (or linearized) dynamics, 
we now seek the corresponding response diagram in our 
highly nonlinear chain. 

Results. The central result is summarized in Fig. 
This represents the “nonlinear spectrum” of the system 
in response to its periodic actuation from the left bound¬ 
ary. To obtain this diagram, a combination of a fixed 
point algorithm (Newton’s method for boundary value 
problems in time) and pseudoarclength continuation was 
used to numerically compute the branches of period-A: 
solutions for various integers k. Similar to m, this was 
accomplished by identifying fixed points of the strobo¬ 
scopic map with a period of k2'K/uj (i.e., k times the 
period of the applied harmonic excitation). The inset 
captures solely the stable (observable) response of the 
system, in the form of the maximal force exerted on the 
wall by the last bead over the period of the solution, 
plotted as a function of the drive frequency. It bears five 
clearly distinguishable peaks (which we will call “reso¬ 
nances” in the sense of local maxima of the response of 
the granular crystal), indicating that at these frequen¬ 
cies a maximal force is transmitted to the right end of 
the crystal. The diagram also features dips (which we 
will call “antiresonances”) between the peaks, associated 
with frequencies that locally minimize the transmitted 
force. This resonance structure is clearly reminiscent of 
a linear spectrum, despite the fact that the basic ingre¬ 
dients of the latter are absent. Moreover, the frequency 
range where these resonances and antiresonances occur 
lies within the nonlinear pass band of the chain HE], i.e., 
the frequency range where solitary pulses or spatially ex¬ 
tended waves can propagate through the corresponding 
granular crystal of infinite extent, whereas for frequencies 
above 2 kHz, there is negligible force transmission. As 
we will see below, the effective phase variation (defined 
below and in the figure captions) between the applied 
excitation and the nonlinear waves that result from it 
and travel along the chain is crucial in the formation of 
these persistent resonance dynamics. We then delve into 
a more detailed exploration of these features. 
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FIG. 1: (Color Online) Maximum force exerted on the fixed wall by the last bead over one period of the periodic motion of 
the crystal. Branches corresponding to periods that are different integer multiples of the period of the actuator are shown in 
different colors. In the quasiperiodic case, the motion becomes dense on a torus, and the average force is estimated over a 
long time frame. The main graph showcases the different branches of nonlinear solutions (denoted by solid when stable, and 
dashed when unstable), while the inset incorporates solely the stable, observable response, i.e., the “nonlinear spectrum”, that 
would be experimentally measured. The amplitude of the periodic excitation used here is 5 x 10~^ m. 


We seek exact periodic orbits for the different frequen¬ 
cies within the range of our actuator (here, between 0 
and 2.5 kHz). The resulting limit cycle solutions (given 
that this is a driven, damped system) are shown in the 
bifurcation diagram of the main part of Fig. There, it 
can be discerned that there is a “backbone” of harmonic 
response (consisting of period-1 nonlinear solutions, with 
oscillation frequencies matching the driving frequency). 
In addition to identifying the existence of such a branch, 
we also explore its stability by means of Floquet analysis, 
identifying the corresponding Floquet multipliers of each 
point; see, e.g., mn] for a detailed discussion. When 
these multipliers are within the unit circle, they corre¬ 
spond to stable dynamics, while their exiting the unit 
circle is associated with instability, which, in turn, is de¬ 
noted by dashed lines in Fig. [T] 

While the period-1 solution is stable at low and high 
frequencies, indicating a robust dynamical response of 
the system in these ranges, it undergoes a series of desta¬ 
bilizing bifurcations in intermediate frequency regimes, 
which lead to the observation of stable branches of higher 
period (subharmonic) solutions, i.e., of periodic responses 
with periods that are integer multiples of the period of 
the applied excitation. These subharmonics are denoted 


by different colors in Fig. We note the existence of 
a resonant peak “within” each branch of periodic solu¬ 
tions, corresponding to a local maximum of the magni¬ 
tude of the largest (over a single period of the solution) 
transmitted force to the right end of the chain. On the 
contrary, antiresonances corresponding to local minima 
of the largest transmitted force typically occur on sta¬ 
ble parts of the period-1 solution branch. The detailed 
spatio-temporal nature of the periodic responses of the 
chain in resonances and antiresonances is of particular 
interest and will be discussed below. From Fig. [l] we 
deduce that the emergence of the period-2 solution oc¬ 
curs through a period-doubling bifurcation (and its dis¬ 
appearance through the “inverse” of the same bifurca¬ 
tion). However, the rest of the branches are isolated (so- 
called “isolas”) and are not emerging directly from the 
primary 1:1 branch. It is also observed, however, that 
there are regions where we did not find low-period sta¬ 
ble periodic solutions; in these regions (see the brown 
dash-dotted lines in Fig. [^, we expect the system re¬ 
sponse to be quasi-periodic (i.e., consisting of two har¬ 
monics with irrationally related frequencies) - or periodic 
with much higher integer periods, as these are known to 
arise in a systematic, interspersed pattern in quasiperi- 
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odic regimes [IMH- 

Figure depicts representative spatio-temporal varia¬ 
tions of the kinetic energy in the regime of the 1:1 and 
1:3 resonance peaks, together with the corresponding dis¬ 
placement waveforms for every bead during one period. 
Visual inspection of the kinetic energy space-time plots 
clearly suggests pulses traveling along the chain (part of 
the “nonlinear spectrum building blocks” mentioned in 
the abstract). These apparent waves are initiated on the 
left end of the chain, travel to the right, become (par¬ 
tially) reflected there, come back to the left end of the 
chain, and there interact again with the next “hit” of the 
actuator. The third component of each figure summa¬ 
rizes this pulse traveling in the form of an effective phase 
variation in the response of each bead: what is plotted is 
the time difference between the first occurrence during a 
period (at stationarity), of positive velocity for each bead, 
and for the leftmost bead of the chain; the almost per¬ 
fect straight line plot clearly implies the traveling pulse 
nature of the long-term, stationary dynamics in the reso¬ 
nance peaks considered. This provides some insight into 
the nature of the relevant periodic, nonlinear solutions 
that dominate the response at resonance of the crystal 
lattice under the periodic external excitation. 

In the parameter regime of the resonance peaks, it ap¬ 
pears that the effective phase variation between the ap¬ 
plied excitation and the traveling pulse initiated in the 
crystal is such that strong energy transfer from the ex¬ 
citing source occurs upon impact. This can be rational¬ 
ized when we take into account that all resonances occur 
within the nonlinear pass band of the crystal, where trav¬ 
eling pulses constitute the nonlinear mechanism of energy 
transmission, and the applied excitation is in the form of 
a periodic train of impulses EHilll]. In the case of 1:1 
resonance (Fig. [^a)) a strong impulse excitation occurs 
at precisely the time instant when the traveling pulse re¬ 
flected from the right immovable boundary reaches the 
left end; one can deduce that the traveling pulses fully 
synchronize with the excitation source. At higher-order 
resonances, however, the time needed by the traveling 
pulse to fully traverse twice the length of the crystal is an 
integer multiple of the period of the applied pulses; for ex¬ 
ample in the 1:3 resonance peak (Fig.j^b)) three periods 
of the excitation correspond to one period of the traveling 
pulse in the crystal; so the excitation inflicts a “strong im¬ 
pulse” at every third period (whereas secondary, weaker 
impulses also occur in between, but they are not strong 
enough to initiate new traveling pulses in the crystal). 
We conclude that the l:n resonance peaks correspond to 
traveling pulses in the crystal with periods equaling n pe¬ 
riods of the applied excitation. This is confirmed in the 
plots of the displacement waveforms of all beads during 
one full period (Fig. [^c) and Fig. [^d)). The traveling 
pulse character of these resonant responses of the crystal 
is also demonstrated in Fig. dfe) and Fig.[|;f), where the 
effective phase variation (as defined above) of the indi¬ 


vidual bead responses is depicted as a function of space 
(more precisely, of bead number). 

Having explored the resonances of our nonlinear spec¬ 
tral picture, let us now turn our attention to the corre¬ 
sponding antiresonances, corresponding to local minima 
of the (largest over one period) transmitted force to the 
right, immovable boundary. As mentioned previously all 
antiresonances are associated with 1:1 periodic responses, 
and this is confirmed by the results depicted in Fig. 
where the second and third antiresonances are shown. 
What is interesting is that these minima of the largest 
transmitted force reside in the pass band of the crystal 
and alternate between resonances. This strongly suggests 
that the observation of antiresonances is not due to the 
incapacity of the granular crystal to transmit energy (in 
the form of traveling solitary pulses) but rather it is due 
to negative interference effects between left- and right¬ 
going traveling pulses; in such cases one might expect to 
observe nonlinear standing waves. This is confirmed by 
the results of Fig. where negative interference effects 
between apparent traveling pulses can be observed, giv¬ 
ing rise to peculiar periodic responses. Indeed, it appears 
that there is a form of co-existence of traveling pulses (as 
in the resonance cases discussed previously) and apparent 
nonlinear standing waves in different parts of the crystal 
(traveling on the left side, standing on the right side, the 
side of the immovable wall). Such states of co-existence 
of different waveforms have been dubbed chimeras in 
other contexts; see for a relevant discussion. This 
co-existence of apparently traveling and standing waves 
highlights the essentially nonlinear nature and the high 
complexity of the long-term, stationary response, and is 
corroborated by the results presented in Fig. [^e) and 
Fig. if). The negative interference of traveling pulses 
propagating in opposite directions in the crystal would 
explain the minimization of transmitted force in antires¬ 
onances, and the “low-intensity” impulses delivered by 
the exciting source in these cases. Such negative inter¬ 
ference appears to lead to the formation of standing wave 
breathing patterns, such as the ones evident on the right 
side of the lattice in Fig. ie) and Fig. |^f). 

All of the above analysis has been performed for a fixed 
value of the forcing amplitude. However, our system per¬ 
mits exploration of a wide parameter space by means of 
the concurrent variation of both the forcing amplitude A 
and the forcing frequency u. This more global perspec¬ 
tive of the system is provided in Fig. Here, we use the 
same color as before (in Fig. to denote the bifurcation 
loci of the different branches. For instance, the isolas of 
the period-4 and period-5 solutions terminate in pairs of 
saddle-node bifurcations (cf. Fig. that are followed, 
in a two-parameter bifurcation diagram as the ampli¬ 
tude A varies, to produce resonance horns within which 
such solutions are contained in the {A, oj) space. We also 
highlight the locus of Neimark-Sacker period-1 bifurca¬ 
tions, giving rise to quasi-periodic solutions. The tips of 
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FIG. 2: (Color Online) (a),(b) Long-term, stationary spatio-temporal variation of the kinetic energy (scaledx 10^) for (a) 1:1 
resonance (340 Hz) and (b) 1:3 resonance (1130 Hz); the input harmonic displacement and applied force excitations for each 
case are also shown on the left of each plot. (c),(d) Displacement waveforms of every bead during one full period of the crystal 
motion. (e),(f) Effective phase variation of the individual bead responses (computed by the time instants when each bead’s 
velocity first becomes positive and when the leftmost bead’s velocity first becomes positive); this highlights the “traveling 
pulse” nature of the long-term, stationary dynamics of the granular chain at the resonance peaks. 


these resonance horns lie, in fact, on the Hopf (Neimark- 
Sacker) curve, when the Floquet multipliers cross the 
unit circle at 4th and 5th roots of unity respectively. 
Needless to say that a host of other resonance horns 
(corresponding to every rational number) and many sec¬ 
ondary bifurcations, both local and global, that cannot 
be feasibly computed, can also, in principle, be found in 
such a diagram (e.g. [MS]). The bottom panels pro¬ 
vide a sense of the variation of one of the system Poincare 
sections (here, by monitoring the plane of the displace¬ 
ment and velocity of the 6th bead). The existence of 
a stable period-5 solution in both panels is illustrated 
by the green points; a coexisting unstable (saddle-type) 
period-5 is shown by the red points, while the unstable 
period-1 solution of the “backbone” branch is shown by 
the isolated point in the middle. The saddle- and node- 
period-5 solutions were initially born in a saddle-node on 
the torus; a remnant of this can be seen in that alternate 
sides of the unstable manifold of the saddles are attracted 
to different stable period-5 points. Global bifurcations, 
involving homoclinic manifold crossings and the associ¬ 
ated horseshoe dynamics (see e.g. [H]), also occur within 
these resonance horns; while in the right panel a coex¬ 
isting secondary stable quasiperiodic attractor emerges 


in the form of an invariant circle. In this last panel, 
one side of the saddle unstable manifolds is attracted by 
the stable period-5, while the other side is attracted by 
this torus; this change clearly suggests that a global bi¬ 
furcation, namely a homoclinic interaction of the saddle 
invariant manifolds with each other, must occur along 
the path from point A to point B in parameter space. 
This is but a glimpse of the considerable complexity of 
this highly nonlinear system and the role of the nonlin¬ 
ear building blocks (stable and unstable periodic, as well 
as quasi-periodic solutions) naturally arising in the non¬ 
linear analogue of the spectral picture we have outlined 
here. 

Conclusions and Future Challenges. In summary, in 
the present work, we focused on forced resonance spectra 
of highly nonlinear homogeneous granular crystals in the 
complete absence of linear acoustics (sonic vacua). De¬ 
spite the lack of sound propagation, in the sense of clas¬ 
sical acoustics, in the medium, we observed that a purely 
nonlinear analogue of the spectrum arises. This consists 
of harmonic and subharmonic responses to the external 
boundary drive of the granular crystal, leading to a se¬ 
ries of peaks (resonances), reminiscent of the transmis¬ 
sion frequency peaks in a finite, linearly coupled lattice. 
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FIG. 3: (Color Online) (a),(b) Long-term, stationary spatio-temporal variation of the kinetic energy (scaledx 10^) for (a) the 
second antiresonance (860 Hz) and (b) the third antiresonance (1,260 Hz); the input harmonic displacement and applied force 
excitations for each case are also shown on the left of each plot. (c),(d) Displacement waveforms of every bead during one full 
period of the crystal motion. (e),(f) Effective phase variation of the individual bead responses (computed by the time instants 
when each bead reaches maximum positive displacement and when the leftmost bead does); this highlights the chimera-like, 
partially “traveling pulse” and partially “standing wave” nature of the long-term, stationary dynamics of the granular chain at 
antiresonances. 


Here, the building blocks of the long-term, stationary re¬ 
sponses are periodic and quasi-periodic orbits, whose ex¬ 
istence and stability we explored computationally. More¬ 
over, we illustrated how the periodic orbits corresponding 
to the energy transmission (resonant) peaks emerge out 
of periodic oscillations within the chain that can be at¬ 
tributed to highly nonlinear traveling pulses; the dips be¬ 
tween them (antiresonances) are due to the emergence of 
a mix of apparent traveling pulses and standing waves, 
trapping the driver’s energy in the vicinity of the lat¬ 
tice edges (a surface breathing mode, resembling a highly 
nonlinear analogue of the well-known Tamm states [25]). 
These results were extended into multiparametric varia¬ 
tions (amplitude as well as frequency of the drive) and 
their generic nature was conhrmed. 

This proof-of-principle demonstration paves the way 
for systematic future investigations of such systems. On 
the one hand, it would be interesting and relevant to at¬ 
tempt to obtain a detailed experimental blueprint of the 
full nonlinear spectrum; indeed, experiments guided by 
these predictions and confirming their validity are un¬ 
derway, and will be reported elsewhere. On the other 
hand, numerous questions deserve further exploration. 


regarding the specifics of corresponding spectral pictures 
in systems with different types of nonlinear couplings, 
or in the presence of local nonlinearities as, e.g., in the 
cradle model of jUj, as well as in the case of internal 
resonators m- A systematic understanding of the mod¬ 
ifications of the spectral characteristics as N progres¬ 
sively increases would be another important aspect for 
future study. Extending these considerations to multi¬ 
dimensional systems with different geometric character¬ 
istics dSl [5S| would also represent a significant general¬ 
ization of the present work. Such studies are presently in 
progress and will be reported in future publications. 
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FIG. 4: (Color Online) (a) Partial two-parameter bifnrcation diagram with respect to the forcing amplitnde and the freqnency 
of the actuator. The curves denote different bifurcations involving the nonlinear periodic solutions. More specihcally: (1) 
the continuation of a Neimark-Sacker bifurcation of the period-1 solution (deep blue); and (2) continuation of saddle-node 
bifurcations of the period-4 and period-5 solutions (cyan and purple, respectively). The saddle-node bifurcations enclose two 
resonance horns that are born on the Neimark-Sacker curve when the Floquet multipliers cross the unit circle at the 4th and 5th 
roots of unity respectively, (b) Zoom-in of the period-5 resonance horn including a one-parameter cut across it. (c) Poincare 
cut showing that the unstable manifolds of the saddles at point A above, along with the stable period-5 solutions, form a closed 
invariant set; at the horn boundary this was an invariant circle on which the saddles and the stable period-5s were “born” 
locked. The axes here are the dimensionless position and velocity of the 6th bead, (d) Moving towards point B in parameter 
space above, homoclinic crossings of the unstable and stable manifolds of the saddles have occurred and two attractors now 
coexist, a stable invariant circle (a quasi-periodic solution) and a stable node period-5 solution. 
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Appendix. We consider a one-dimensional chain of 
spherical elastic beads in Hertzian contact with no pre¬ 
compression: 


TO- 


d^ui 

dt"^ 


TO- 


TO- 


d^Ui 

dt'^ 


d?UN 

dt^ 


i^*\/ ¥ (^TT/t) - ■Ui)+'^^ - (mi - M2)+^ 

+D [(27r/A' cos {2Trft) — Ui)H{A' sin (27r/t) — ui) — (ui — U 2 )H(ui — U 2 )] , 

4^* n3/2 , s3/2 

3^ \ ~2 “ '“i+i)+ 

+D [(ui_i - Ui)H{u^_i - Ui) - [iii - Ui+i)H[ui - Ui+i)] 

(wAf-i - un )+‘^ - { un )+'^ 

+D [(uiv-i - un)H{un_i - un) - {un)H{un)] , 


r 


where 

TO 

R 

A' 

f 

D 

E* 

E 

V 

Define 

Xi 


w = 27r/, 


— mass of each bead (28.84 g) 

— radius of each bead (9.525 mm) 

— amplitude of the actuator (5 x 10“^ m) 

— frequency of the actuator 

— viscous damping coefficient 
between beads (100 Ns/m) 

E 


A = — 


2(1-zz2) 

Young’s modulus (193 GPa) 
Poisson’s ratio (0.3). 


A! (original amplitude of the actuator) 

Uj 

^0 
A' 

A-o 
rt 
D 


where 








3to 


Using the variables defined above we end up with our 
nondimensionalized system of equations 


Xi = {x^-i - xffl'^ - {Xi - 

-LA [(ij-i - Xi)H{xi_i - Xi) 

with xn^i = xn-\-i = 0 and xq = Asin{ujt). 
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